COMPUTING BORCHERDS PRODUCTS 
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Abstract. We present an algorithm for computing Borcherds products, which 
has polynomial runtime. This improves on the currently common method, 
which has exponential runtime. An implementation of the new algorithm 
shows that it is superior also in practice. 



1. Introduction 

Product expansions are a utile construction in the world of automorphic forms. 
For example, writing r and r' for two variables in the Poincare upper half plane, 
we can express the famous j-invariant as follows: 

]{T)-j(T) = {e -e ) [I (l-e 1 ') 

m,n=l 

for certain integral coefficients c mn . This formula holds for sufficiently large imagi- 
nary parts of t and t'. Several facts can be deduced from it. For instance, one sees 
immediately that close to infinity, the j-function is one-to-one on the fundamental 
domain. 

It turns out that product expansions exist for a much larger class of automorphic 
forms, namely orthogonal modular forms for 02, n (R). Thanks to the Koecher 
principle, such product expansions allow for a direct investigation of all zeros of 
such an automorphic form, regardless of the limited domain of convergence they 
have. 

The convergence of these more general product expansions remained an open 
problem for a long time. When Borcherds presented his first work on product ex- 
pansions for orthogonal groups |Bor98| (see also |Bor95| IHM981 IKon97| ) . he started 
a revolution. Product expansions, which at that time played a role in the classi- 
cal theory of elliptic modular forms and Jacobi forms only, became important to 
other areas in mathematics. One branch of research aimed at Kac-Moody algebras 
which could be constructed by means of lattice theory [Gri95, GN98, Sch04, Sch08]. 
As a natural generalization of finite dimensional Lie algebras, they attracted the 
physicists' attention. On the automorphic forms side, research in this area focuses 
on orthogonal groups for rather big lattices. A second, equally popular branch of 
research originating in Borcherds' work is concerned with the structure of divisor 
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class groups and Chern classes |Bor99| iBorOOl IBru02| [BOTOl ICG08j . It is founded 
on the fact that one can easily read off the divisor of an automorphic form given by 
a product expansion. Lattices with particular nice structure, which are typically of 
moderate size, are in the focus of mathematicians working in this direction. The 
Fourier expansion of Borcherds products is a particular aspect that is interesting 
in this context. Among other things, it can be used to understand the structure of 
graded rings of modular forms |DK04| IKri05|, |MaylO| , or to check for relations like 
the ones coming from Saito-Kurokawa-Maafi-Gritsenko lifts [HM1Q]. Most impor- 
tantly, the multiplicative version of the degeneracy of BPS dyons, occurring in string 
theory, are Borcherds product of this type (see, for example, |HM98llJS05[|DM llj). 
The last named author will return to this subject in the sequel. 

Large scale computational approaches to Siegel modular forms were pioneered by 
Skoruppa in |Sko92]. His considerations became famous, in particular, because in 
the course of his studies, he discovered two exceptional Hecke splittings in weight 24 
and 26, which until today remain unexplained. Equally important computations 
that led to deep insights have been performed by Ibukiyama |Ibu02| . and Poor 
and Yuen |PY00[ [QPY08| . A recent, striking result was obtained in |PY09| . where 
strong evidence for the Paramodular Conjecture [Yos80l rYos07] . which parallels the 
Shimura-Taniyama Conjecture, was provided using paramodular forms, a variation 
of Siegel modular forms. 

The purpose of this work is to provide an algorithm for explicit computations 
with Borcherds products that runs in polynomial time. Such an algorithm is vastly 
faster than the direct evaluation of Borcherds' formula, which requires repeated 
multiplication of increasingly large polynomials, leading to an exponential runtime. 
The presented algorithm is based on two ideas: to make the problem as linear as 
possible and to split off parts that behave distinctly. In order to achieve the first 
objective, we apply expo log to Borcherds' original formula, and we evaluate the 
logarithmic term first. This avoids the repeated multiplication of polynomials. To 
achieve the second objective, we analyze the structure of Fourier indices that are 
positive with respect to a fixed Weyl chamber (see (|3.4p to (|3.8|l ). This allows us 
to reduce the number of terms taken into consideration. 

We demonstrate the actual superiority of our approach by providing an imple- 
mentation for Hcrmitian modular forms over Q(V=3) in Sage |S + llbl iDGPSlO] . 
We compare it to a usual implementation: for moderate precision (see Section [4] for 
a discussion of precisions) the new algorithm is already faster by a factor between 2 
and 4, depending on the memory architecture of the used hardware. The coefficient 
of the leading term of the new implementation's runtime is seemingly large. This 
prevents it from improving even more on the runtime of the naive implementation. 
It originates, however, in the not optimal memory management of Sage, and it will 
reduce drastically for a reimplementation in C or CH — h 

In Section [21 we review the basic theory of orthogonal modular forms and 
Borcherds products. Section [3] deals with the algorithm to compute Fourier ex- 
pansions of Borcherds products. The implementation for Hermitian modular forms 
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is discussed in Section |4l In particular, this section contains a comparison of run- 
time tests. 

2. Preliminaries 

2.1. Orthogonal modular forms. Throughout the paper, we fix an even lattice 
L of signature (2,n) with n > 3. The quadratic form attached to L is denoted by 
q, while we write (/x, A) for the bilinear form q(/j, + A) - q(fJ.) - <z(A). If A is a vector 
in a lattice with fixed and obvious embedding into L<S>Q, we understand that q(X) 
is the quadratic length of A therein. The dual lattice of L is 

L* = {A e L ® Q : (fi, A) e Z for all fi e £,}. 

Clearly, L£ L*. Thus we can construct the Z- module L&/L, that is equipped with 
the quadratic form q taking values in Q/Z. The pair (L#/L, q) is a finite quadratic 
module of order \L\. It is called the discriminant module of L, and we denote it by 
discL. 

We write O(L) for the orthogonal group associated to L. It is the group of all 
Z-linear bijections of L that leave the quadratic form q invariant. Any such trans- 
formation acts also on and discL. The stable orthogonal group O(L) [disci] 
consists of all transformations fixing discL elementwise. 

We consider a fixed connected component T> + e of 

V e := {A € L <8> C : q(X) = A (A, A) > 0}. 

The projectification T> + = P(2?+), which is a Hermitian locally symmetric domain 
of type IV, is the natural domain of definition for orthogonal modular forms. We 
define them using the -k- homogeneous pullback to T> e . A function / : T> e -* C is 
called -fc-homogeneous if /(sA) = s- k }{\) for all s e C and all A e D e . 

Definition 2.1. An orthogonal modular form of weight k for a finite index subgroup 
r < O(L) is a holomorphic, -k -homogeneous function f : T> e -> C satisfying the 
following conditions: 

(1) For all 7 e T we have f °J = f, 

(2) For every isotropic vector A e ihoCD + e the function f(y\) is bounded as y -*■ oo. 

After fixing a section T> + ^ T> ei we freely identify any orthogonal modular form 
with its restriction to T> + . 

The Borcherds products that we are interested in are arithmetic orthogonal mod- 
ular forms. Throughout the paper, we assume that L = U ffi U © (-l)Lo, where U 
is the unimodular hyperbolic lattice of rank 2 and Lq is a fixed even positive def- 
inite lattice. Notice that discL = disc(-l)Lo- Under this assumption, there is a 
canonical section T> + ■-»■ T> e whose image has typical elements (l, q(z) - tt' , t, t', z), 
t, t' e C, z e Lq <S> C. In this setting, any orthogonal modular form has a Fourier 
expansion 

V a([a, b, c]) e(aT + cr' + (b, z)), 

where e(x) := e 2mx . The sum ranges over triples [a, b, c] with a, b e Q and 6 e L <S>Q. 
We write Qq for the additive group of such triples. The submonoid of integral 
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indices [a,b, c] € Qq that satisfy a,c e Z and b e L*, is denoted by Q. All Fourier 
expansions that we will deal with satisfy a([a,6,c]) = 0, if [a, 6, c] £ Q. 

We call an index positive definite or semi-definite if disc([a, 6, c]) := ac-q(b) > 
and a > 0, or disc([a, 6, c]) > and a, c > 0. We write [a, 6,c] > if [a, 6, c] is 
positive definite, and [a, b, c] > if it is positive semi-definite. The monoids of 
integral positive semi-definite indices is denoted by Q + . 

The group M := 0(U © (-1)L ) c 0(U ®U® (-1)L ) gives rise to symmetries 
of the Fourier expansion. That is, we have a([a,6,c]) = x( m ) cx(m[a,b,c]) for all 
me M and a fixed character x of M. 

Remark 2.2. For the time being, drop the condition L ~ U © U © (— l)£o ; 
we /laue imposed on L. If n > 5, it is always possible to find a sublattice L' £ L 
of finite index that splits: L' a uiC/ ffi © (-I)Lq wzi/i positive constants u\ 
and U2 (see, e.g., [Bor98, Section 8]^. Since the orthogonal groups of L and L' 
are commensurable, the considerations in this paper apply to arbitrary Borcherds 
products, if only n > 5, even though the details are tedious to work out. 

2.2. Borcherds products. We revise the construction of Borcherds products, 
which dates back to |Bor95) and |Bor98) . A more accessible discussion, valid in 
some interesting special cases, can be found in |Gri95| IGN 98J. Given an elliptic 
vector-valued modular form, we construct an orthogonal modular form. 

Let H := {t e C : Jm(r) > 0} denote the Poincare upper half-plane. Since we 
deal with half-integral weight, the next definition involves the metaplectic cover 
Mp 2 (Z) of SL2(Z). It is the preimage of SL2(Z) in Mp 2 (R), the connected double 
cover of SL 2 (R). The primages of ker(SL 2 (Z) -» SL 2 (Z/iVZ)) in Mp 2 (Z), where 
N runs through N, will be called principal congruence subgroups. 

Write g = { a c b d ) for a typical element of SL 2 (R). The elements of Mp 2 (R) can 
be written (g, r h> \/ct + d), where the first component is an element of SL2(M) 
and the second is a holomorphic function on H. Since there are two branches of 
the square root, this yields indeed a double cover of SLa(M). 

Given a representation (p,V p ) of Mp 2 (Z), k e |Z and F : H -> V p , we define 

(F\ kp (g, Vc-Td)) (r) :- JZFTd* p((g, V^d)) f(^). 

Definition 2.3. Let (p,V p ) be a finite- dimensional representation o/Mp 2 (Z). A 
weakly holomorphic vector-valued modular form of type p and weight k e ^Z is a 
holomorphic function F : H -*■ V p such that the following conditions are satisfied: 

(1) For all g e Mp 2 (Z) we have F\k. p g = F. 

(2) The Fourier expansion of F has the form 

Y, /(n)e(nr), f(n) e V p . 

-oo«neQ 

Remark 2.4. If the weight k is integral and p is trivial on some principal con- 
gruence subgroup, the representation p and the factor of automorphy factor through 
SL 2 (Z). 

The most relevant case to us is p = p^, where p^ is the Weyl representation 
associated to the finite quadratic module discL. The space V p has a canonical 
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basis index by discL. We refer the reader to (Bor98, Sko08] for a definition and 
more details. Suppose that F is an elliptic modular form of type pl- The Fourier 
expansion of its components F\ is 

Fx(r)= £ /(A,»)e(nr) (2.1) 

-00«TI 

with A e discL. By /(A,n) with A e L# we mean f(X,n). 

With a fixed basis w±,..., w n -2 of Lq <g> K given, we write b > if there is 1 < 
j < n - 2 such that (b,Wj) > and (b,Wj') = for all f < j. A triple [a,6,c] is 
called positive, [a, b, c] > 0, if c > 0, or c = and a > 0, or a = c = and 6 > 0. This 
definition is motivated by Weyl chambers discussed in |Bor981 Section 6]. We can 
find a Weyl chamber W for U ffi (-I)Lq such that ([a, 6, c],W) > if and only if 
[a, 6, c] > 0, where the Wj are rational vectors. 

Borcherds also defined a Weyl vector W F attached to F and W. In our setting 
(see also |Der01| ). we have: 

tt w := ^ E /fr-dis^)), bw '' = \ ^ /( fo ' _disc ( 6 )) 

OO 

== «w - E "^™) E /(&,-»- disc(6)), 

: = e(affT + c^r' + (bw, zyj. 

With this notation at hand, we reformulate [Bor98, Theorem 13.3]. Note that 
we suppress a root of unity in the definition of ff. 

Theorem 2.5. Fix a weakly-holomorphic vector-valued elliptic modular form F 
of type pl and weight (2 - n)/2 with Fourier expansion (12. ip . Then the following 
product converges and is an orthogonal modular form for 0(L)[discL]: 

V F :=W F [I (l- e (ar + cr' + (6,z))) /(Misc([QAc])) . (2.2) 

[a,b,c]>0 

Crucial to our approach is the function 

$F=l g* F = ] Dg W F + E f(b,[a,b,c])log(l-e(aT + CT'+(b,z))), (2.3) 

[a,b,c]>0 

which is well-defined and holomorphic in the Weyl chamber W. We have = 
exp($ F ). 

3. An algorithm for Fourier expansions of Borcherds products 

We write Q[Q] for the direct product of 1-dimensional Q- vector spaces indexed 
by Q. Any element will be written as a possibly infinite sum of multiples of e^ a ' b ' c ^ 
for [a,b,c] e Q. 

Fix a vector- valued elliptic modular form F of type pi, and weight (2 - n)/2. 
Define 

^ oo n[a,6.c] 

®f ■■= - E f( b > disc (K 6 > c])) E — - 6 Q[Q1- (3-i) 

[a,b,c]>0 n=l n 
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Clearly, <&f >-»■ $f - log Wf as e^ a ' b ' c ^ is mapped to e(ar + ct' + (b,z)). 

Even though Q[Q] carries no algebra structure, we can define powers of 
Write (j) F ([a,b,c\) for the coefficient of e ^ b ' c \ We set 

£ ( E n^([^.^^])) e[a ' M » (3-2) 

[a,b,c] Y}i.i[ai-bi,Ci] = [a.b,c\ 

where the [a^, &i, c^J's are elements of Q. By the next lemma, the inner sum is finite. 
Lemma 3.1. Given I e N, the set 

|([a i ,& l ,c J ]). =1 n e Q l : = [a,b,c], (^([a*, &i, c;]) * /or i} 

z=l 

is finite. 

Proof. For reasons of symmetry, it suffices to prove that there are only finitely 
many [ai,i>i,ci] that possibly contribute, regardless of the [a,,0j,Cj]'s with i + 1. 
Assume that the tuple ([ctiTbijCi])^ e Q l is an element of the above set. 

Since F is a weakly holomorphic modular form, there is a lower bound d such 
that f(b,D) = 0, whenever D < d. By definition of the relation > and of <£>/, 
only indices [ai,6i,ci] satisfying c\ > contribute. In particular, we can find 
an upper bound on c\ for those [oi,&i,Ci] that contribute. Consider those with 
ci > 0. If <Pf([cli, 0i, ci]) + 0, we can write [oi,6i,ci] = m[3i,6i,ci] with m e N 
and disc([2i,6i,ci]) > d. This gives rise to the inequality disc([ai, bi, Ci]) > c 2 d. 
From this, we deduce lower and upper bounds A\ and ^4 U on a\ for those [oi, b\, c{\ 
with ci > that contribute. For every fixed a%, the same inequality excludes all but 
finitely many b\, since Lq is definite. 

In the case of c\ = 0, only indices [oi, b±, 0] satisfying a\ > contribute. By what 
we observed in the case c\ > 0, we have a\ < -(I - l)min{0, A,} + a. That is, we 
have an upper bound on a\ for all [ai,6i,ci] that contribute. 

Assume that a\ > 0. Using the factorization [oi,6i,Q] = m[5i,6i,Q] analogous 
to the one above, we conclude that there are only finitely many &i's that can occur. 

Consider the case a\ = c\ = 0. Then bi > 0. Recall the rational vectors 
Wx, ... , w n -2 6 Lq ® Q discussed in Section [2] Choose j between 1 and n - 2 such 
that (wj,bi) > and (vjj',bi) = for all 1 < j' < j. Since Wj is rational, the set 
of possible (wj,bi) for b\ > is discrete in K. Because Y,i(wj,bi) = (vjj,b) is fixed 
and, in particular, independent of b\, there is only a finite number of &i's that can 
occur by the same argument as above. □ 

The following theorem is used implicitly, when computing Borcherds products, 
whether or not Borcherds' formula is evaluated only naively. Nevertheless, no proof 
seems to be available in the literature. 

Theorem 3.2 (Formal Borcherds convergence theorem). The exponential 

* F := W F exp(<f>F) = W F Y ( 3 -3) 

is well-defined in the sense that, given [a,b,c], the coefficient of e^ a ' b ' c ^ in <fr F is 
nonzero for only finitely many n. 
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Moreover, for any [a,b,c] the coefficient of e^ a:b ' c ^ in $^ equals the Fourier 
coefficient of e(ar + cr' + (6, z)) in 



Proof. The first part follows by analyzing the proof of Lemma I3TT1 in the light of the 
new statement. The second part follows from the first part, since by Theorem l2.5[ 
due to Borcherds, the product expansion for &p converges locally uniformly abso- 
lutely. □ 

We split $f into five pieces: 



2t:= £ -f(b,disc([a,b,c])) £ — — , (3.4) 

[a,b,c]>0:a,c>0 n=l 



<8:= £ -f(b,disc([a,b,c]))Z—^> ( 3 - 5 ) 

[a,b,c]fO:a>0,c>0 



Y -f(b,dwc([a,b,c]))Z—r- ' ( 3 - 6 ) 

[a,b,c] : a<0, c>0 



©:= Y, -f(b,disc([a,b,c]))Y,—^> ( 3 - 7 ) 

[a,b,c] : a>0, c=0 



E -/(&,dfec([a,6,c]))£ ■ (3- 

[a,b,c] :a=c=0,b>0 n=l n 





[a,6,c] 




n 


e" 


[a,f>,c] 




n 


e n[ 


a,b.c] 




n 


e n[ 


a,6,c] 




n 




a,6,c] 





We can further split € = Y.]=i £f 



e 



<* r .= Y -f(b,di S c([a,b,c]))Y 

[a,b,c] :a=c=0,b>0 n=l " 

where the sum only ranges over those b satisfying (wj,b) > and (wj',b) = for all 
1 < f < j- Note that 21 is M-invariant, while the four other pieces are not. 

We use capital letters to denote elements of Q[Q] and lower case letters to denote 
their coefficients. For example, we write £([a, 6, c]) for the coefficient of e^- a ' b ' c ^ in 

In the last step of the next algorithm, the multinomial coefficients 

i + - + nA (Sj=inj)! 



/ni + ••• + m\ _ { Lj=i 
\ nx,...,ni ) ft 1 



for m, . . . , ni 6 N will occur. 

Algorithm 3.3. The following algorithm computes the Fourier expansion $>p for 
alle [a ' b ' c] with a,c< B e N: 
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1: Set d <- max{d : f{b, d') = for all d' < d and all b}. 

2: Set Oneg < — (B - cw ~ 1) d. 

3: For < f3 < min{F + a ncg - aw - a, B - cw - ct} 

and < 7 < B - cw - a do Step 4 to 6 and then jump to Step 13. 
4: Compute <- ^ £ 7 . 
5: Set a m - ln «- min{a : 3 [a, 6, c] : £([a, 6, c]) 0}. 
6: For all 5 < B - a m i n do S'tep 7 and 8. 

7: Set S^ j7 .,5 <- 3,3. 7 J)" 5 and 5 <- S^^j f/or abbreviation only). 

8: For 1 < j < n - 2 recursively do Step 9 to Step 12 and, at any leaf of 
the associated decision tree, set 3^ 7 <- 3 wi£/i r\ = (ry 1; . . . ,77^-2)- 

9: Set b m in,j «- min{(iyy, 6) : 3[a,6,c] :£([a,6,c]) #0}. 
if: iSef & TO ax,j ■ max{(wj, 6) : 3 [a, 6, c],£([a,6,c]) : [a, 6, c] > 0}. 
11: For all 0<T]j < (-b min ,j + 6 maxj ) / min{(wj , b) > : 6 > 0} do Step 13. 
12: Set 5 ^ 3 <*f . 

13: Compute £ (a + + 7 + <* + J?)!" 1 2t Q H ft7ltfj1J . 

ITie /ast sum ranges over < a < B - a w and all 3^ ~ g v that were defined. 

Proof of correctness. The correctness follows immediately, when making the bounds 
used in the proof of Lemma 13-11 explicit . □ 

Remark 3.4. Considering each step separately, one easily sees that Alaorithm \3.3\ 
has polynomial runtime in B. Evaluating formula (|2.2j) directly has exponential 
runtime, since the coefficients of a weakly holomorphic modular forms grow expo- 
nentially. 

4. BORCHERDS PRODUCTS FOR THE HERMITIAN MODULAR GROUP 

The authors provide an implementation of Algorithm [33] for Hermitian modular 
forms over the imaginary quadratic fields Q(-n/D). It is written in Sage [S + llb| . 
and the authors think of it as a proof of concept. It is equally based on another 
implementation presented in [Raul lb] by the last named author, providing a model 
for general Fourier expansions. This dependency results in restricted functionality 
for D + -3; For all these cases the framework currently provides only limited sup- 
port. Both implementations are available on the authors' homepage |Raulla| and 
are intended for prompt integration into Purple Sage [S + lla| . 

As a basis for the maximal order O £ Q(\AD) we use 1 / \AD, (1 + VD) / 2, 
writing b = bi / yD + b^il + ^/D) / 2. With an appropriate Weyl chamber chosen, the 
condition b > 0, described in Section |2~2| is equivalent to 02 < v (62 = a b\ < 0); 
A condition easy to test. When computing the powers of 21, the implementation 
makes use of its GL2(0)-invariance. Crucial to the implementation is the fact 
that is the Fourier expansion of a Hermitian modular form. This allows for a 
significant reduction of the number of coefficients that need to be calculated. For 
more details on the implementation, the reader it referred to comments to be found 
in the source code. 
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TABLE 1. Time needed to compute the coefficients of 045 for indices 
with diagonal entries bounded by B using a stub implementation. 



B 


5 


6 


7 


8 


9 


Naive evaluation 


0.04 s 


7.15 s 


502.24 s 


12349.70 s 


43.5 h 


New algorithm 


0.53 s 


15.48 s 


271.00 s 


5029.15 s 


31 h 



To demonstrate the value of our approach, we compare our implementation with 
the common, naive, method. The example that the authors computed is ^45, the 
weight 45 Hcrmitian modular form for Q(v-3) (see, e.g., |Der01| ). For the sake 
of clearness, the authors have implemented the algorithm in Sage. The memory 
footprint of Sage, however, is larger by a factor exceeding 40 compared to native im- 
plementations, and the memory access is not optimal. For this reason the provided 
implementation suffers from cache misses and an enormous memory consumption 
for higher precisions. In Table [TJ the impact of Sage's memory consumption can be 
observed. For precision 8 and 9 the performance gain realized up to then increases 
much slower than expected. The reason is that for example in the case of preci- 
sion 9, the memory management of Sage (presumably due to non-efficient iteration 
over large lists) allocates 100 GB of memory, while theoretically the new algorithm 
does not consume more memory than the naive evaluation does. Nonetheless, the 
runtimes for precisions 7 to 9 show that the new algorithm is not only asymptot- 
ically faster, but for moderate precisions yields a clear advantage over the naive 
one. Note that the examples in Table [T] were computed using a Sun server of the 
Lehrstuhl A fur Mathematik at the RWTH Aachen University. The relative in- 
crease in performance was even stronger when using other systems, based on AMD 
chipsets, with more efficient memory architecture. For lack of memory these could 
not be used to run all tests. 

We conclude that an implementation of the presented algorithm aiming at real 
applications must be implemented in C, C++ or a language alike. This will be a 
comparatively easy task with the provided Sage code at hand. 
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